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ABSTRACT 

Utilizing -four-quadrant experimental data determined by Minio- 
vich converted to the modi-fied advance and an g 1 e-o-f -advance 
coe-f-f icient -form, complete Karman-Knapp Circle Diagrams are 
plotted -for a number o-f 3 bladed propellers. Reducing the 
-family o-f curves to a rational basis, algebraic equations are 
derived to model the per-formance o-f a three bladed propeller 
through all -four quadrants using Eulerian turbomachine theory. 
The result is a set o-f equations which can be used to generate 
complete propeller thrust and torque characteristics given only 
incomplete per-formance data suitable -for use in steady state 
design or dynamic simulation. 
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CHAPTER 1 



INTRODUCTION 



1 . 1 Primary motivation -for this kind of model ino method 

For' a naval architect or anyone interested in the modeling o-f 
propellers as a part o-f the design procedure leading to an 
integrated ship design, the problem long persisted where one 
is required to estimate the performance of a propeller for 
ship modeling purposes before the design has progressed to the 
point where one is reasonably confident as to what the final 
desired propeller should look like. Even more frustrating is 
the problem of having somehow determined the geometry of the 
propel 1 er but not having sufficient experimental data to allow 
a comprehensive computer simulation. 

Even where it has existed, considerable effort has been expen- 
ded to convert the data into useable form for computer simu- 
lations since the practice of tabulating the results as func- 
tions of an advance coefficient has difficulty in the vicinity 
of zero ship or shaft speed. 

The problems associated with modeling a machine such as a 
propeller are formidable since a complete model must make 
allowances for complicated fluid motions which are only just 
beginning to be understood. In particular, it is still quite 
difficult to pass completely from propeller geometry to the 
detailed lift and drag forces required for specific propeller 
char act eristics. 

However, there now exists a considerable amount of documen- 
tation for modeling machines very much like a propeller. For 
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example, industrial axial turbines and pumps have had the 
attention oi modelers and those who perform simulations -for 
decades. Taken together, propellers, turbines and pumps are in 
■fact, members o-f a general class o-f machines known as turbo- 
machines. Models -for many o-f the large industrial machines 
have been quite success-ful in predicting machine per-formance 
characteristics and are widely used in control system synth- 
esis and design. A general model -for propellers which could 
incorporate that data base would have a -firmer basis in exper- 
ience and -first principles. 

1 .2 Some current propeller modeling methods 

At present there are two common methods used by designers to 
obtain propeller data -for computer simulations. The -first 
involves experimental determination o-f the characteristics o-f 
a physical model o-f the -full size propeller -followed by con- 
version o-f the data into the -familar thrust and torque coe-f-fi- 
cient -form. Such coe-f-f icients are then used in a table look-up 
scheme within the simulation program to obtain the necessary 
thrust and torque values. There are advantages to this kind o-f 
method. First, the determination o-f the propeller characteris- 
tics through experimentation (scaling considerations aside) 
ensure that the data used in the simulation will be as accu- 
rate as possible. Secondly, the use o-f a model o-f the intended 
propeller acts as a test o-f the geometry and reveals at an 
early date any serious conceptual problems. 

The disadvantages o-f such a method are -fairly obvious. To test 
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a model o-f a propeller at an early stage in a design, you must 
have a well de-fined idea as to the range o-f operati on o-f the 
propeller and its desired geometry. To allow -for a comprehen- 
sive simulation, you must test and record the model data in 
all -four quadrants o-f operation. Since cost is always a signi- 
-ficant consideration, the machining o-f a propeller model and 
the set-up and operational costs o-f the tests represent a 
sizeable investment in the design. As alluded to earlier, raw 
data is not suitable -for computer simulations due to the 
problem o-f 'singularities' o-f the advance coe-f-f icient at zero 
ship speed and zero sha-ft rpm. The experimentally determined 
data must be converted to a -form that avoids such di-f-ficul- 
ties. 

The second method o-f modeling propeller characteristics invol- 
ves complex and relatively recent work in propeller hydrody- 
namics. The computer models o-f Pro-fessors Kerwin C13 and 
Cummings [2] are quite thorough and provide some revealing 
results. As a 'complete' model, these computer programs are 
undoubtedly at the vangard in propeller theory. With a -firm 
basis in hydrodynamic theory, the continual evolution o-f these 
programs will ultimately achieve very accurate results. Un-for- 
tunately, these programs are immensely long and require consi- 
derable amounts o-f CPU time. Because o-f the attendant expense, 
they are not yet practical -for incorporation directly into a 
ship simulation program. 
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1 .3 The Objective o-f the Modeling Equations 



The objective of this thesis will be to derive a set of rela- 
tively simple algebraic equations that will serve as a tool 
for naval architects in the ship design process. By being able 
to predict accurately propel 1 er performance in all four quad- 
rants of operation based on Knowledge obtained from limited 
experiments and in a form suitable for computer simulation, 
the naval architect can save considerable time and expense in 
the simulation portion of the design. Experimental determina- 
tion of the points of minimum dissipation would be sufficient 
to allow for a comprehensive simulation for forward or crash- 
back maneuvers. Such a set of equations may not eliminate the 
need for actual experiment in the final stages of the design, 
however they would greatly simplify the initial requirements. 
Additionally, where experimental data exists for particular 
propellers in only the first quadrant, the equations could be 
used to reconstruct the remaining three quadrants with good 
accuracy . 
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CHAPTER 2 



DERI»v¥i»TIQN o-f the ALGEBRAIC MODEL 



2. 1 Fundamentals o-f Eulerian turbomachines 

The practice o-f presenting open water propeller per-formance 
data as nondimensi onal thrust <Kt) and torque <Kq) coe-f-fi- 
cients versus an absolute advance coe-f-f icient (lambda) or 
simply advance coe-f-f icient (J) C33,C43 incorporates assump- 
tions which inherently restrict the machines o-f interest to 
Eulerian turbomachines. That is, by considering only high 
Reynolds numbers the per-formance o-f the machine is governed by 
inertia and pressure -forces. Dimensional analysis on the ex- 
pression : 

(Thrust) F=-f (Vp,D,N,p) 



yiel ds 



Kt= F/pD-»N2=-f(^.»p/ND) =-f(J) 



( 2 . 1 ) 



Similarly -for torque, 

Kq= T/PD5N2=g((.>p/ND) =g(J) (2.2) 

where Vp is the speed o-f advance o-f the propeller through the 
water and D any characteristic propeller dimension (usually 
diameter). This single parameter representation o-f propeller 
per-formance yields the power-ful result that the per-formance 
data o-f geometrically similar propellers will, when plotted 
together, generate to good approximation a single curve. 

The same assumptions would , -from a -fluid mechanics stand- 
point, give Euler^s equation o-f inviscid motion. That is. 



pvXGRAD(v) = -GRAD(P) 



(2.3) 
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Or, in a more general -form to incorporate unsteady -flows, 



PD<v)/Dt = -6RAD<P) <2. 4) 

where D/Dt is the substantive derivative. 

In conjunction with equation (2.4), we express mass conserva- 
tion by; 

D(p)/Dt + PGRAD3(v =0 (2.5) 

where P is independent o-f position and v is the vector -fluid 
velocity. Re-ference [53 chapter 3 and Re-ference C63 are use-ful 
for a more detailed analysis. 

A machine satisfying equations (2.4) and (2.5) will, as veri- 
fied exper imen tal 1 y , demonstrate what is referred to as 
■Eulerian similitude". That is, geometr i cal 1 y similar machines 
will produce similar performance characteristics when present- 
ed in a properly chosen nondimensi onal form. Casanady [73 

refers to this as a 'performance law'. 

Based on the discussion above, the determination of the per- 
formance char ac ter i st i cs of a large turbomachine is usually 
done with small scale models. The full scale results are 
derived from the model results by recognizing that at corres- 
ponding points on a charac ter i st i c curve, all nondimensi onal 
parameter ccmbi nat i ons are also equal. 

However, there are problems with scale effects not evident in 



such a single parameter representation . In addition to geo- 
metric dissimilarities, such as roughness and clearances, hy- 
drodynamic di-f-ferences are also present. These include Rey- 
nold'^s and Froude number variations as well as cavitation and 
compressibility number di-f-ferences. Attempting to incorporate 
some or all of these variations into a performance law would 
require that the law be multidimensional. As an example, one 
could attempt to formulate a nondimensi onal representation for 
thrust and torque characteristics by: 

Kt=f(J,Re) Kq=f<J,Re) . 

Here, the equivalence of the advance coefficient in conjun- 
ction with the Reynold'^s number would lead to similitude. 
Similarly, one could also form yet another expression for 
thrust and torque 

Kt=f<J,P> Kq=f(J,P> 

where P represents a nondimensi onal ratio for propeller pitch. 
Such a formulation could be applicable to Kaplan turbines or 
controllable pitch propellers where the blade pitch is varied 
mechanical 1 y . 

Multi-dimensional performance laws unfortunately lead to fami- 
lies of curves, a problem which single parameter laws avoid. 
As an example, Bovet ESI presents efficiency versus output at 
a constant head for a Kaplan turbine for various blade set- 
tings. An analytical treatment for a multi-dimensional per- 
formance law in terms of generating numerical data would 
entail a matrix of possible states with dimension equal to the 



12 



number o-f nondimensional characteristic parameters. Obviously, 
dimensional analysis techniques used to generate characteris- 
tic nondimensional per-formance parameters are limited in their 
applicability in per-formance simulation. To allow incorpora- 
tion o-f -fluid or geometrical a-f-fects into the analysis would 
require another approach. That approach should produce a for- 
mulation which retains the advantages of similitude yet allows 
for some fluid or geometrical effects. 

2.2 Ideal, canonical and perfect ship propeller equations 

Figure 2.1 shows a propeller blade cross section at some 
arbitrary radial distance from the propeller hub under the 

condition of zero nominal speed of advance. Figure 2.1a in- 

dicates the tangential velocity of the blade element <ki<*r> and 
fluid velocity <Vf) . Figure 2.1b shows the lift <L) and drag 
(D) forces with respect to an effective angle of attack a. The 
magnitude of the effective angle of attack is a function of 
the circulation which causes an induced velocity £31. 

If we define a torque exerted by the element as Fr then, 

Fr= LSin<cc) * DCos<o:) <2. 6) 

Further, defining Ff as a force perpendicu 1 ar to Fr in the 
direction of fluid motion (fluid pressure term), 

Ff= LCos(a) - DSin(tx) (2.7) 
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<a) 




<b) 



Figure 2 . 1 
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Since L and D are expressed as: 



L= Cl<J4PAvXv) ; D= Cd(J^PAvXv) 



<2. 6) and <2. 7) can be rewritten; 



Fr= ^4PAv<C13CV-f + Cd3«Vr) <2.8> 

F-f= ^4PAv<Cd3fVr - Cl XV-f) (2.9) 



In per unit (see Appendix A), (2.8) and (2.9) become: 



m= v(ClXq + CdXn) 
p= v(ClXn - CdXq) 



Or using the canonical notation, 

m= Rg3(q + RdXn 
p= RgXn - RdKq 



( 2 . 10 ) 

( 2 . 1 1 ) 



( 2 . 12 ) 

(2. 13) 



where 



Rg=rqXv 

Rd=rdXv 



rq=rq(P) v=SQRT(nXn+qXq) 

rd=rd(P) ^ARCTAN(q/n) 



The correspondence with the bond graph representati on shown in 
Appendix B is now readily apparent. 

The mechanical power can be expressed as FrXVr while the -fluid 
power can be expressed as F-fXV-f. Thus, dissipation losses can 
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be -found to be 



Pd= FrJfVr - Ff5^V-f 



<2. 14) 



In per unit, (2.14) becomes 



Pd= mXn - pXq 



<2. 15) 



Substituting (2.12) and (2.13) yields 



Pd= RdXvXv 



(2. 16) 



Hence, as would be expected, an ideal ship propeller would 
have no power lost in dissipation i -f it could generate li-ft 
without drag. The equations reduce to the rather trivial 
expressi ons 



with Rg representing an energy trans-f ormat i on modulus. Rd can 
be considered as a general loss term that includes a combin- 
ation o-f a number o-f loss mechanisms such as -form drag and 
separation losses. Writing Rd in a more general -form typical 
o-f such losses such that at the point o-f minimum dissipation 
loss, Rd is identically zero 193, we can de-fine a PERFECT 
propel 1 er . 



m = Rg^q 



(2. 17a) 



p = RgSfn 



(2. 17b) 
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Rd = W(^)X<<n-q)3f<n-q)/v) 



(2. 18) 



Using the Kutta condition to -formulate a general expression 
•for the modulus Rg would suggest 



Rg = rg X< ^) X < <nX n) -< qX q) )/v <2.19) 



di-f-fering -from the -form proposed by Paynter 
which could represent some linear function in 
<2.18) and <2.19) into <2.12) and <2.13) and 
arbitrariness in W<P) and X<P) , we can arrive 
expressi ons: 



by the term rg 
n. Substituting 
recognizing the 
at the f ol 1 owi ng 



m = rgXq * K10)Xv<n-q) 
p = rgXn * K2<F)Xv<n-q) 



< 2 . 20 ) 

< 2 . 21 ) 



where K1 and K2 now incorporate the individual character of a 
given propeller. The equations above for a PERFECT propeller 
<zero minimum dissipation loss) can be au^ented to represent 
a PRACTICAL propel 1 er by including an additional term to 
reflect the actual non-zero minimum dissipation loss. 



m = rgXq * K10)Xv<n-q) 
p = rgXn K2<P)Xv<n-q) 



+ RlXn 
- R2Xq 



< 2 . 22 ) 

<2.23) 



R1 and R2 would be functions of the advance coefficient J for 
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a given propeller, but in a simpler model could be reduced to 
constants. As will be seen later, the values o-f R1 and R2 have 
a direct e-f-fect on the shape o-f the 4-quadrant curves o-f a 
given machine and in particular, the slope o-f the zero torque 
and thrust radial s. K10> and K2(P) -for an idealized propeller 
would be identical. The degree o-f agreement between Ki and K2 
based on experimental data would be a measure o-f the validity 
o-f the model -for a given rg. 

2.3 Determination o-f the algebraic modeling equations 

Using equations (2.22) and (2.23) as a model -for determining 
propeller characteristics requires a determination o-f the K(^) 
term as well as Rl, R2 and rg based on experimental data. The 
data produced by Miniovich CIO] and converted to a modi-fied 
advance coe-f-f icient -form Cll] was used in an abbreviated 
methodology in order to -formulate a general model that could 
be expanded upon in a more detailed analysis. See also Re-fer- 
ence C 12] . Such an analysis would include an investigation o-f 
those -factors which collectively yield the individual charac- 
ter o-f a given machine yet have de-fineable speci-fic in-fluen- 
ces, the details o-f which are beyond the scope o-f this thesis. 
In Re-ference Cll], experimental results -for thrust and torque 
are tabulated -for eighteen three bl aded propellers with vary- 
ing expanded area ratios and pi tch-to-di ameter ratios. The 
propeller designations are those used by Baker and Patterson. 
Table 2.1 is provided as a quick re-ference -for EAR and P/D 
values -for each propeller. 
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Figures B3 and B4 in Appendix B are plots o-f raw data -for two 
series o-f propellers. To be use-ful , these curves must be 
rationalized to a consistent basis. To do so, the points o-f 
minimum dissipation were -found -for each o-f the propellers (See 
Appendix A> . Those points were then used to scale -figures B3 
and B4 such that the point o-f minimum dissipation in the 
region o-f normal pump operation (quadrant I> was -fixed at 
( 1 , i> and in the region o-f normal turbine operation (quadrant 
III) at (-1,-1) in the q vs n plane. The mapping o-f the 
data to those points o-f minimum dissipation is done over the 
entire q vs n plane recognizing the ■'mirror image' pattern o-f 
the curves. That is, the point o-f minimum dissipation in the 
-first quadrant -for positive torque and thrust is re-flected 
into the third quadrant on the negative torque and thrust 
curves. Hence the scaling is re-flected through the origin 
rather than simply across the zero torque and thrust radials. 
Figures 2.2 and 2.3 are plots o-f the thrust and torque curves 
respectively -for three propellers selected -for their varying 
expanded area and pi tch— to-di ameter ratios and are representa- 
tive o-f the -form o-f the 4-quadrant curves -for all o-f the 
propellers. With the data plotted in this -form, the similarity 
o-f the torque and thrust curves versus the advance angle is 
striking. In addition, the general -form o-f the curves is quite 
consistent in the regularity o-f the speci-fic perturbations 
which in-fluence all o-f the propellers in the same manner. 
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Figure 2.2 
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Figure 2.3 
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ProD 


No . 


1 


2 


3 


4 


5 


6 




EAR 


0.5 


0.5 


0.5 


0.5 


0.5 


0.5 




P/D 


0.6 


0.8 


1.0 


1.2 


1.4 


1.6 


ProD 


No . 


7 


8 


9 


10 


11 


12 




EAR 


0.8 


0.8 


0.8 


0.8 


0.8 


0.8 




P/D 


0.6 


0.8 


1.0 


1.2 


1.4 


1.6 


ProD 


No . 


13 


14 


15 


16 


17 


18 




EAR 


1.1 


1.1 


1. 1 


1 . 1 


1 . 1 


1. 1 




P/D 


0.6 


0.8 


1 .0 


1.2 


1.4 


1.6 



Table 2. 1 

Before continuing, I must point out that there is an 
approximation which was used which requires an explanation. As 
Figure B2 in Appendix B demonstrates, the lines of zero 
thrust (head) in the first and second quadrants are obviously 
not 180 degrees apart. The difference between the angle 
measured clockwise from the first quadrant zero thrust radial 
and the second minus 180 degrees I call the skew angle. The 
skew angle for high specific speed axial flow pumps is 
generally quite small. See reference [133. The skew angle 
computed for the propellers used in this analysis were all 
less than 3 degrees and were generally less than 1 degree. The 
amount of skew angle is a function of the magnitude of losses 
of the machine. Since the skew angle is quite small, the 
curves were approximated to a PERFECT propeller by dividing 
the skew angle between the two radial s found from the minimum 
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dissipation calculations. The scaling values for q and n were 
then recalculated using these adjusted angles. As a result, 
the first and third quadrant asymptotes were 180 degrees 
apart. The approximation was felt Justified since the calcu- 
lated minimum dissipation angles were obtained with a confi- 
dence of plus or minus 2 degrees due to the required extrapol- 
ations in the correspondi ng portion of the experimental data. 
Such an approximation also simplifies the derivation of the 
modeling equations considerably. Table 2.2 lists the 
calculated minimum dissipation angles and the approximate 
angles in both quadrants for six of the eighteen propellers 
along with the calculated and approximated magnitudes of the 
thrust and torque radials at those angles. The approximate 
values are shown in parenthesis. 

Prop. No. Min. Dissipation Radial Length 





angl es 


thrust 


torque 


1 


35 . 2 < 34.6) 

-146.0 <-145. 4) 


1 1 .52< 8. 97) 
3.92<3.80) 


15.39< 14.88) 
25.83<22.44) 


2 


39.5 < 38.5) 

-142.5 <-141.5) 


5.64<5.02) 

4.98<4.59) 


13.54< 11.32) 
25.65<23.88) 


3 


44.5 < 44.3) 

-136.0 <-135. 8) 


4.95<4.82) 

5.02<4.90) 


10.54< 9.96) 
19.99< 18.68) 


5 


53.5 < 53.0) 

-127.5 <-127. 0) 


5.06<4.78) 

6.88<6.31) 


9.27< 8.85) 
26.99<20.92) 


15 


45.6 < 45.6) 

-134.5 <-134. 5) 


6.99<6.90) 

4.82<4.79) 


12. 10< 12.04) 
20.36< 19.61) 


18 


58.0 < 57.8) 

-122.5 <-122. 3) 


6.99<6.59) 

7.48<7.07) 


11.96< 11.47) 
20.99< 18.76) 



Table 2.2 



Using equations (2.20) and <2.21> for the PERFECT propeller 
and setting rg equal to n, we can plot the K(P) functions 
required to fit the experimental data by setting m and p to a 
constant (1) and solving for K1 and K2 at each data point. 
Since K1 and K2 must be positive at all times, we account for 
the sign change in m and p in the different quadrants as the 
calculations proceed. By plotting K1 and -K2 as functions of 
the normalized advance angle, we can gain a quick qualitative 
feel for the similarity of K1 and K2. Figures 2.4, 2.5, and 
2.6 are plots of K1 and -K2 for propellers 9, 15 and 18 
respec t ivel y . 

As is readily apparent, even with a simplified function for rg 
<n) , K1 and K2 are very similar. The behavior of the K<P) 
functions in response to changes in pi tch-to-di ameter and 
expanded area ratios is also apparent and are manifest in the 
amplitude variations and relative smoothness of the plots. 

On the basis of that similarity, we form a new pair of K(^) 
functions by the strai gh tf orward transformation; 

Ka(P) = (Kl + K25/2 

Kb(P) = (Kl - K2)/2 

such that Kl equals Ka+Kb and K2 equals Ka-Kb. Hence, Ka 
represents an average of the two functions while Kb is a 
measure of their difference. Figure 2.7 is a plot of the Ka 
and Kb functions for propel 1 er number 15. 

A different approach may be to account for the general ampli- 
tude differences in Kl and K2 on eitherside of the zero thrust 
and torque radial s and themselves. Four gain parameters could 
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Figure 




26 



Figure 
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Figure 




<a) 




Figure 2.7 
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be obtained -for each propeller such that when Ka and Kb are 
plotted, the maximum amplitude of KbO> will be considerably 
reduced. This was done for propeller number 15 by gain scaling 
the average amplitudes as outlined above. Figure 2.8 shows the 
result for Ka and Kb. KbO> now approaches what may be con- 
sidered as individual propeller nuance characteristics. 

Since KaO) is essentially doubly periodic over the range of ^ 
we form yet another function that can be used to approximate 
K1 and K2 valid over the entire range of ^ . By overlaying the 
plots of KaO) for -135 < P < 45 and 45 < P < 225 and 
again averaging we obtain an approximate K(^> function applic- 
able to both the thrust and torque curves and valid for the 
two ranges of P stated above. 

To retain all of the information available such that we are 
able to reconstruct the original K1 and K2 functions as 
necessary, we must define four additional functions. Those 
functions are formed by first separating Ka and Kb on the 
intervals of -135 < P < 45 degrees and 45 < ^ < 225 degrees. 

Labi ed: 



Kal = Ka 
Ka2 = Ka 
Kbl = Kb 
Kb2 = Kb 



-135 < P < 45 degrees 
45 < ^ < 225 degrees 
-135 < P < 45 degrees 
45 < P < 225 degrees. 



We then define the four functions that represent the average 
and differences of these functions on the two intervals of P. 
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(Kal+Ka2)/2 



K = 

1 

K = <Kal-Ka2)/2 
2 

K = (Kbl+Kb2)/2 
3 

K = <Kbl-Kb2)/2 



4 



Algebraic manipulation will show that: 

K1 = K +K +K +K -135 < P < 45 deg 
12 3 4 

K1 = K -K +K -K 45 < P < 225 deg 

12 3 4 

K2 = K +K -K -K -135 < P < 45 deg 
12 3 4 

K2 = K -K -K +K 45 < P < 225 deg. 

12 3 4 

Our approximate K(P) -function is now -found -from: 

K<P) = (K -i-K )/2 
1 3 

valid on the intervals of -135<^45 and 45<^225 degrees. 

If the two functions K1 and K2 were the same and the identical 
underlying function was responsible for the thrust and torque 
characteristics on each side of the 4-quadrant asymptote, then 
K , K , and K would be zero. In summary, K(P> is an average of 



1 3 4 

K1 and K2 which is then symmetrized for the -•• and - torque and 
thrust curves. Figure 2.9 shows the result for propeller 
number 9 on the interval -PI < P < PI . 

To determine some of the underlying fundamental characteris- 
tics of the K<P> function we can formulate an expansion of 
K<P) . Two of the methods which are available involve the Walsh 
functions and the Fourier series. With the Walsh functions, 
square waves are superimposed and in the limit can be used to 
approximate simple sinusoids. They could be particularly use- 
ful here since the K(F) function does have very sharp minimums 
at the 45 and -135 degree angles. The Fourier series is s«ne- 
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( a) 




(b> 



Figure 2.8 
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Figure 



what more -familar and was used here. As a dual to the Walsh 



■functions, in the limit, the Fouries series can be used to 
approximate square waves. 

By obtaining the Fourier series expansion -for K<P) we have 
•finally, a potential model -for this -family o-f propellers. 
Table 2.3 lists the -first -five Fourier coe-f-f icients -for 
propellers 9, 15 and 18. The Fourier expansion is taken over 

the interval between the minimum dissipation angles. At those 
angles the K<P> function approaches zero and provides a 
convenient point to begin the expansion. The expansion is 
given by equation (2.24> where c is defined by equation 
<2.25) . 

By substituting the K(P) expression into the basic modeling 
equations (2.20) and <2.21> , we can compare the experimental 1 y 
determined thrust and torque curves to the model by plotting 

the equations over the range of -1 < n < •►! . In doing so, we 

are assuming that the differences in K1 and K2 are sufficient- 
ly small to permit their omission. Justification of this 
assumption lies in the accuracy of the curves generated by the 

model equations when compared to actual experimental data. 

Figures 2.10 and 2.11 shcxM the results for propeller number 
18. 



KO) = aO SUM< a<k)Cos<kx) b<k)Sin<kx) ) <2.24) 



c = 



SQRT< aJ^a + b3<b ) 



<2.25) 




Figure 2. 10 
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Figure 2.11 



35 



Propeller Number 







9 






15 






18 






a 


b 


c 


a 


b 


c 


a 


b 


c 


o 

II 


5.20 


0 


5.20 


6.27 


0 


6.27 


10.31 


0 


10.31 


k=l 


.59 


-.26 


.64 


.64 


.59 


.87 


-.23 


2.98 


2.99 


k=2 


-2.02 


- 1.06 


2.28 


- 2.63 


-1.15 


2.92 


-3.27 


-1.29 


3.52 


k=3 


-1.25 


.27 


1.28 


-.57 


-. 19 


.60 


-1.25 


. 12 


1.26 


II 


-.59 


-.87 


1 .06 


.78 


-.53 


.54 


-.78 


.01 


.78 



Table 2.3 



To demonstrate the power o-f the model the curve resulting -from 
utilizing only the zeroth term o-f the Fourier series expansion 
is plotted along with the curve plotted using the -first three 
terms. Doing so reveals that the di-f-ficulty in -fitting the 
proper K -function lies essentially in the dissipation regions 
o-f the 4-quadrant curves. In a sense, this is expected since 
the KO> -function represents a modulating dissipation term in 
the basic modeling equations. What is interesting is that the 
plots indicate that a common dissipation mechanism a-f-fects 
both the torque and thrust characteristics and are lumped 
together in the KO) term. 

Additionally, it is obvious that quite accurate results -for 
normal ahead propeller operation and -for turbine operation can 
be achieved with KO) assumed to be a constant. As more terms 
o-f the expansion are added, the -fit can be improved in the 
dissipation regions. For a loss o-f power simulation -for exam- 
ple, the most accurate results -for thrust and torque values 
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would require a large number of terms. The limitation on the 
number of terms that can properly be added -to the expansion 
lies in the residue of the K1 and K2 functions. Recalling that 
K(F> was formed from an averaging process involving K1 and K2, 
the functions which represent their differences were formed to 
allow retrieval of the original K1 and K2 data. To the extent 
that these differences degrade the KO) approximation, their 
impact should properly be less than the last term added to the 
expansion. Again , the Nalsh functions may reduce the number of 
terms required particularly in the normal pump and turbine 
areas of the 4-quadrant diagram. The model equations can also 
be used to construct the familar Kt and Kq versus J curves. Kt 
is given by equation <2.26) and Kq by equation <2. 27) . J is 
simply q/n where R= n3<n+qXq and THETA=Arctan J. 

2 

Kt = l./RXcos (THETA) (2.26) 

2 

Kq = l./RXcos (THETA) (2.27) 

Here, the R and n refer to the appropriate values for the 
thrust and torque curves. 

By selecting an n and solving the modeling equations (2.20) 

and (2.21) for q with a given K(P) , Kt and Kq can be quickly 

determined for any value of J. Figure 2.12 is a plot of the Kt 
and Kq curves for propeller number 18 using the zeroth term of 
the K<^> expansion and the first three terms for J > zero. 
Note that the scale factors used in constructing the rationa- 
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lized plots o-f thrust and torque must be considered in calcu- 
lating Kt and Kq so that their ranges o-f values are those 
normally encountered in the literature. The scale -factors used 
-for propeller number 18 in constructing Firgure 2.12 can be 
deduced -from Table 2.2. 
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<b) 

Fi gure 2.12 
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CHAPTER 3 



INTRODUCTION o-f ADDITIONAL TERMS and OTHER 
DETAILS in PROPELLER CHARACTERISTICS 



3. 1 Pitch as a primary determinant 

As will be discussed in Chapter 4, the variation in pitch 
expressed as a pi tch-to-di ameter ratio <P/D) has a significant 
effect on the performance characteristics of a given propel- 
ler. From a designer's standpoint, Saunders [33 addresses 
pi tch-to-di ameter selection primarily on the basis of obtain- 
ing maximum propeller efficiency, particularly when utilizing 
series chart data in the initial stages of the design. From 
the perspective of the modeling equations however, it is more 
convenient to see the relationship between the relative pitch 
to diameter ratio and specific speed. References [143 and [153 
present descriptions on the use of the specific speed para- 
meter . 

Specific speed variations have a significant effect on the 4- 
quadrant diagram not easily seen in the propellers presented 
here. Since propellers are in general very high specific speed 
machines, the impact of variations of that parameter is some- 
what obscure. The principle effect of specific speed variat- 
ions on the 4-quadrant diagram was only touched on in the 
plots of the thrust and torque charac ter i st i cs discussed ear- 
lier. It is seen primarily in the location of the zero thrust 
and torque radials of Table 2.2. As the relative specific 
speed increases, the radials tend to rotate counterclockwise 
while the skew angle addressed in Chapter 2 tends to decrease. 
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How variations in pitch a-f-fect the specific characteristics of 
a propeller is best seen in an analysis of the performance of 
controllable, reversible pitch propellers <CRP> . Reference 
t 163 is quite thorough in examining quantitatively the impact 
of pitch on the efficiency and thrust and torque coefficients 
of a CRP propeller. It is seen in the modeling equations as a 
dominant factor in the determination of the amplitude of the 
KO> function. An increase in the pi tch-to-di ameter ratio has 
two noticeable effects. First, it tends to increase the rela- 
tive amplitudes of the various harmonic components in K<^) 
which results in a somewhat wider fluctuation of the 4- 
quadrant curve in the principle dissipation regions. Secondly, 
it has a direct affect on the zeroth term of the Fourier 
expansion which represents the average value of the KO) fun- 
ction. The net effect is a shrinking of the 4-quadrant curves 
(corresponding to larger torque and thrust coefficients). 
Obviously, the 'dc' value of the expansion is dominant since 
an examination of Figures 2.2 and 2.3 would indicate that the 
harmonic content of K<P) in the dissipation regions, specifi- 
cally in quadrant IV is apparently proportionally equal for 
all three propellers. 

3.2 Number of blades and blade area 

The effect of variations in numbers of blades was not examined 
in detail. However, as pointed out in Chapter 4, the formula- 
tion of the modeling equations was quite general and utilizes 
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what is recognized as the actuator disk theory to obtain their 
initial relationships. A comparison o-f 4-quadrant diagrams -for 
three b1 aded and seven bladed propellers has not indicated any 
potential problems. Nhat is really the primary determinant is 
the amount o-f blade area expressed here as expanded area 
ratio. Yet Figures B4<a) , B4<b> , 2.2 and 2.3 would seem even 
to discount any major in-fluence o-f EAR on torque and thrust 
charac ter i sti cs in general. 

3.3 E-f-fects o-f cavitation 

Cavitation is probably the most di-f-ficult aspect o-f propeller 
operation with which to deal. Rightly so, a lot o-f time and 
energy is being expended to describe the mechanics o-f cavita- 
tion and its in-fluence on propeller charac ter i sti cs. Re-feren- 
ces £4], £143, £173, £183, and £193 explore qualitatively and 

quantitatively the role o-f cavitation in propeller perform- 
ance. Re-ference £143 relates the speed o-f advance o-f a propel- 
ler to the cavitation number. By de-fining the cavitation 
number <sigma> with equation <3.1), 

sigma=2g<h-^hA-hy)/Va3^'v'a <3.1) 

where g is the acceleration due to gravity, h the depth o-f the 
propeller axis beneath the water sur-face, hA the atmospheric 
head pressure, hV the vapor pressure o-f water and Va the 
propeller speed o-f advance, it is shown that <using typical 
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values) si^na is related to Va by the relationship of equation 
(3.2) . 



sigma= 1070/Va3(Va (3.2) 

with Ma in knots. How the cavitation number affects Kt and Kq 
is seen in Reference C4], Chapter 16, Figure 108 as a fall off 
of Kt and efficiency as the advance ratio decreases. It does 
appear that the effect of changes in sigma for a higher J is 
smaller, also seen in Reference C 163 . 

Yet nowhere is the role of cavitation included directly in the 
modeling equations. Since the intention of the modeling equa- 
tions is to provide a reasonable method of determining propel- 
l er characteristics at an early stage of a ship design say for 
the purposes of simulation, perhaps an extension to include 
cavitation affects is unwarranted. However, we can anticipate 
where such affects may be present. 

Equations (2.22) and (2.23) can be modified to indicate gen- 
eral relationship to cavitation number shown below. 

m=rg(si^a) 3(q + Kl(P,sigma)5(vX(n-q) + Rl(sigma)^n (3.3) 
p=rg(si gma) 3(n + K2( P, si gma) 3(v3( ( n-q) - R2(si^a)^q (3.4) 

In some cases, the presence of cavitation is handled by an 
effective reduction of the propeller blade area. This could be 
incorporated into the K(P) functions. However, the greater 
dissipation would alter the 4-quadrant zero thrust radial s to 
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the R1 and R2 terms are included in this 



some extent, hence, 
case. Notice that R1 is also a function of sigma since the 
losses will be reflected from the hydraulic to mechanial side. 
Although the potential is there for inclusion of cavitation 
affects, such a step may be beyond the scope of the assumed 
model equations. Again, the intent of the equations is to 
provide a relatively simple means to simulate propeller per- 
formance for first-cut design purposes. Any extension beyond 
that goal will increase the complexity of the equations rapid- 
ly. 



3.4 Hu 1 1 /propel 1 er interactions 



Two of the familar parameters which reflect the influence of 
hul 1 /propel 1 er interactions are the thrust deduction factor 
< t) and the wake fraction <w) . The wake fraction appears as a 
consequence of the disparity in ship forward velocity and the 
velocity as seen by the propeller. The wake fraction is de- 
fined as: 



w = <V-Va)/V 

where V is the ship velocity and Va the velocity that the 
propeller sees. 

The thrust deduction factor results from the requirement of 
the propeller to produce enough thrust to propel the ship 
through a resistance of Rt while also reducing the high pres- 
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sure area over the ship's stern. De-fined as: 

t = <T-Rt)/T 

where T is the net required thrust to propel the ship at a 
given speed. Values -for both -factors are usually obtained -from 
model tests. Re-ference C 16] would indicate that t and w vary 

little over a wide range o-f ship speeds. In terms o-f the 

modeling equations, the approach was one o-f considering the 
propeller in an open water condition and hence, corrections 
-for t and w were not included or necessary. Owing to the 

actuator disk approach essentially taken here, allowances -for 
complicated in-flow -fluid geometries due to disturbances ahead 
o-f the propeller have been ignored. 

Re-ference C 16] also addresses the pitch dependency o-f the 
thrust deduction -factor and gives general ranges o-f t -for P/D 
ratios o-f -•■ 1 to -1. Re-ferences C4] and C20] address the incor- 
poration o-f the thrust deduction and wake -fractions in more 
detai 1 . 

3.5 Other considerations 

As -formulated, the modeling equatins simply do not address 

many o-f the -fine details o-f propeller design. Most o-f these 
however, including blade rake, skew-back, the shaping o-f blade 
edges, hub -fairing and blade pro-files are in large measure 
details which go beyond the requirements -for the initial 
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stages o-f propeller design. Once the ''macroscopic' character- 
isics o-f a generic propel 1 er type is determined, the impact o-f 
cavitation and speci-fic blade loading requirements can be 
examined but only in the detail allowed by actual model tests 
or by a complex computer program specializing in propeller 
hydrodynamics. 

A tacit assumption which has been incorporated -fr«n the begin- 
ning is that the nondimensi onal , per unit quantity q is ob- 
tainable -from the dimensional values o-f propel 1 er speed o-f 
advance Va and in turn, the -flow rate through the propeller. 
There is, however, a term called the axial in-flow -factor a [4] 
which is included in the momentum theory -for propeller action 
not directly accounted -for in the calculations. It has been 
assumed to be a quantifiable constant. 

The e-f-ficiency o-f a given propeller and its relationship to 
the behavior o-f the KO) -function has not been explored. The 
calculation o-f the e-f-ficiency o-f a given machine is straight- 
-forward, but the qualitative relationship o-f the e-f-ficiency to 
pi tch-to-di ameter or expanded area ratio in the modeling 
equations is yet to be determined. That the points o-f minimum 
dissipation or points o-f best ef-ficiency can be determined 
-from experimental data or in-ferred -from other considerations 
is accepted, but the use o-f the modeling equations to estab- 
lish a 'best' e-f-ficiency machine given speci-fic operating 
requirements is a topic that begs -further work. 
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CHAPTER 4 



PREDICTION of PROPELLER PERFORMANCE USING THE 
MODELING EQUATIONS 



4 . 1 Comparison with very similar propellers 

Nith only a -few o-f the eighteen propellers available invest- 
igated in detail regarding the specific K<P) functions, never- 
theless it is possible to extend the results to the remainder 
of the family of propellers. As seen by a comparison of pro- 
peller numbers 9, 15, and 18, the K(^) function retains the 
same general shape over the range of P/D ratios of 1.0 to 1.6 
and EAR'S of 0.8 to 1.1. Intuitively at least, the higher EAR 
appears to smooth out the K1 and K2 functions in the dissipa- 
tion regions while slightly increasing the average magnitude 
overall. The effect of EAR is more easily seen in the raw data 
plots of Appendix 6. In Figure B4<b) for example, the higher 
EAR has its most pronounced effect at P equals +90 and -90 
degrees. This is also shown in Figures 2.2, 2.3, 2.4, and 2.5. 
A point by point comparison will also indicate that the varia- 
tion in EAR has little effect on the thrust characteristics 
other than at those two particular points of the 4-quadrant 
curve. Seen more clearly in Figure 64<a> , it is reinforced in 
Figures 2.2 through 2.5. The expanded area ratio has as some- 
what larger effect on the K1 curve (torque), particularly in 
the normal turbine region of the 4-quadrant curve. This is not 
obvious in Figure 84(b) where the plot of the raw data would 
indicate that EAR has essentially no effect outside the prin- 
ciple dissipation regions. 
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The e-f-fect o-f changes in pi tch-to-di ameter ratio is much more 
obvious. Beginning with the raw data o-f Figure B3<a> , we can 
see the e-f-fect o-f P/D ratio throughout the analysis. An exam- 
ination o-f Figures 2.2 and 2.3 show this in-fluence quite well. 
Note however, contrary to intuition, the closer to the origin 
that a given point lies, the higher the relative magnitudes o-f 
the corresponding Ct and Cq coe-f-f icients. A comparison o-f 
Figures 2.5 and 2.6 reveal that the in-fluence o-f P/D ratio is 
nearly the same -for both the thrust and torque characteris- 
tics. Table 2.3 points out that the various principle -frequen- 
cies o-f the K<P) -function change in amplitude as well as the 
average 'dc' term. The substantial variations in the peak 
amplitudes o-f the K<P) -function are also evidenced in the raw 
data o-f Appendix B. The apparent critical points in Figure 
B3<a) are intriguing. The increase in P/D ratio is re-flected 
in the relative speci-fic speeds o-f the various propellers and 
is demonstrated in the gradual counterclockwise rotation o-f 
the zero thrust and torque radial s. Those apparent critical 
points must lie along the pivot axis o-f that rotation but the 
reason -for the correspondence isn't clear. 

By examining Tables 2.2 and 2.3 and Figures 2.2 through 2.6 as 
well as the -figures o-f Appendix B, it is possible to in-fer 
with very reasonable accuracy the appropriate values o-f the 
rating points and the -fundamental values o-f KO) -for all o-f 
the propel 1 ers o-f the data base. Nith that knowledge, the 
thrust and torque plots could be generated and compared with 
Figures 2.2 and 2.3. 
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By plotting the higher harmonic amplitudes versus the zeroth 
harmonic amplitude, a comparison can be made between three of 
the propellers examined in detail regarding the K(P) function. 
Figure 4.1 shows the amplitudes of the first and second 
harmonics plotted against the zeroth harmonic amplitude for 
propellers number 9, 15, and 18. By extrapolating along the 

curves for a propeller whose pi tch-to-di ameter and expanded 
area ratios would suggest a particular zeroth harmonic 
amplitude, we can approximate the values of the Fourier series 
expansion terms for the first two harmonics. 

As an example, propel 1 er number 1 has a pi tch-to-di ameter 
ratio of 0.6 and an expanded area ratio of 0.5. Assuming the 
zeroth harmonic amplitude is equal to four, the first two 
harmonic amplitudes can be inferred by extending the curves to 
the appropriate value of the zeroth harmonic. Using the 
modeling equations with those amplitudes of the respective 
Fourier series expansion terms, a comparison can be made 
against the experimental 1 y determined thrust and torque 
curves. Figure 4.2 shows the result for the positive thrust 
curve . 

Owing to the influence of the higher harmonics near J equals 
zero, some error is expected in that region of the thrust 
curve. However, as Figure 4.2 indicates, even extrapolated 
values of the K<P) expansion terms yield quite reasonable 
results for the normal pump and turbine portions of the 
curves. The value of the zeroth harmonic (first term of the 
Fourier series expansion) has a substantial influence in the 
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Figure 4.2 
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neighborhood o-f J equals zero and acts with the neglected 
higher order terms to produce part o-f the error near the J 

equals zero values o-f Figure 4.2. Ni th more study o-f the 

trends in the harmonic amplitudes with pi tch-to-diameter ratio 
and expanded area ratio, better approximations o-f the thrust 
and torque curves o-f similar propellers will be possible. 
Another important source o-f error in Figure 4.2 is in the plot 
o-f the experimental data. The scaling which was done to 
produce the rationalized data depends on the values o-f n and q 
at the point o-f minimum dissipation. Those scale -factors are 
quite sensitive to errors in the location o-f those points, 
su-f-f icientl y so to cause errors much like those of Figure 4.2. 
A comparison of the curves of Figure 4.2 also points out a 
weakness of the Fourier series representation as opposed to 
other possible expansion methods in forming the KO) approxi- 
mation. The absence of the sharper peaks in the curve gener- 
ated by the modeling equations indicates that a fairly large 

number of terms of a Fourier series expansion is required to 
reproduce the experimental data in the fourth quadrant dissi- 
pation region. A piecewise continuous approximation of K<P) in 
the fourth quadrant would provide improved results for a 
simulation requiring better correlation to experimental data 
during say, transient analysis. 
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4.2 Extrapolation to other ppopelleps 



As developed in this thesis, the modeling equations are quite 
general and although they owe their establishment to a parti- 
cular set of three bladed propellers, they are equally appli- 
cabl e to other types of propellers as well. Hecker [213 pre- 
sents 2-quadrant data for a few seven bladed propellers 
including a hybrid. Plotting the experimental data in the same 
manner as with the three bladed machines, one is struck with 
the similarity of form of the curves to those presented earl- 
ier. Although not shown here, a comparison of the experimental 
data to that generated by the modeling equations using only 
approximate values for the various parameters was quite suc- 
cessful. Although the principle areas of error were in the 
dissipation regions as expected, the comparison did not indi- 
cate any serious problems with the applicabilty of the mode- 
ling equations. A more detailed comparison of the results 
obtained from the modeling equations to much different propel- 
l ers is yet to be done. Considering the scarcity of complete 
4-quadrant experimental data for a large variety of propeller 
types and the limitations inherent in this form of presenta- 
tion, such a comparison wasn^t possible. However, the success 
of the equations to date would tend to enforce the argument 
for their wide ranging applicability. 
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4.3 Extrapolation to missino quadrants 



One of the advantages of the modeling equations lies in their 
ability to reconstruct essentially the entire 4-quadrant diag- 
ram using only partial information. Particularly when suffic- 
ient information is available to determine the point of mini- 
mum dissipation, a "universal" K<P) function can be scaled to 
fit the rationalized data and the appropriate values of ^p and 
ND at the point of minimum dissipation substituted to generate 
either the nondimensional Ct and Cq values or the dimensional 
values of thrust and torque directly. If only the first quad- 
rant point of minimum dissipation is known, the third quadrant 
can be approximated with little loss of accuracy. This is 
particularly useful where first quadrant data exists but the 
model from which the data was taken is no longer available. 
Extending the concept slightly, a comparison of available data 
on the rating points of various propellers could form a data 
base from which the approximate rating point of a specific 
machine could be inferred. Thus, if only very limited experi- 
mental results are available, the 4-quadrant curve could be 
constructed with some confidence by comparing the modeling 
equations results with whatever information is available. 
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CHAPTER 5 



SUMT^RY 



5.1 An overview o-f the method 

The application o-f Eulerian turbomachine theory to the model- 
ing and simulation o-f pump/turbine systems has a long and 
success-ful history. While the particular method used here is 
somewhat more recent, this thesis represents an initial at- 
tempt to apply the methodology to high speci-fic speed extended 
machines. Using very complete experimental data -for a -family 
o-f three bladed propellers with varying expanded area and 
pi tch-to-di ameter ratios as a basis, this procedure was able 
to reveal some of the important fundamental behavior under- 
lying the dissipation mechanisms of propellers in general. 
Perhaps more importantly, the model equations were able to 
point out a common dissipation mechanism for thrust and tor- 
que, i.e., the mechanisms which contribute to the overall 
energy losses of the system can be seen to 'reflect' through 
the transformation of rotational mechanical to translational 
fluid energy. The plotting of the K(^) function for different 
propel 1 ers showed their strong correlation in form, a correla- 
tion which could in a more thorough examination be seen in the 
relative magnitudes of the components of the Fourier series 
expansions. Certainly, the zeroth component of the expansion 
represents an important parameter since a quite reasonable fit 
can be made over much of the four quadrants utilizing it 
solely in the determination of the model 4-quadrant diagram. 
This method also introduces to the simulation of propeller 
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characteristics a procedure which avoids most o-f the 
difficulty with singular values of the advance coefficient. It 
reduces the number of required curves to represent complete 
propeller characteristics to only two, both of which are 
reasonably well behaved and quite smooth. The model equations 
are relatively simple and basically strai ghtf orward. Thus, 
they lend themselves to computer manipulation with little 
di f f i cu 1 ty . 

From the standpoint of the ship designer, the simulation of 
ship dynamics which include a realistic propeller in the 
initial stages of the design would be much easier and certain- 
ly more cost effective. A generic set of equations with the 
specific values of a given propeller rating point incorporated 
in them would be sufficient to provide reasonable results in 
the initial design stages which are oriented more toward near 
steady state ahead and astern operation than a thorough simu- 
lation over all quadrants. Once the geometry of the propeller 
is determined, the equations can be updated to provide in- 
creased accuracy. That accuracy depends on the number of terms 
of the expansion one is willing to incorporate into the simu- 
lation program. Thus, a scale model of the propeller is not 
necessary until the need arises for a self-propelled model 
test or a detailed investigation is required into the cavita- 
tion or vibration characteristics of the propeller. 
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5.2 The potential -for -further applications 



Since this thesis represents an initial examination o-f the 
applicability o-f Eulerian turbomachine theory to propellers, 
little has been addressed regarding a quantitative comparison 
o-f various propeller types. It was demonstrated that the 
results obtained -for the three bladed propellers which -formed 
the data base -for this analysis could be applied to seven 
bladed propellers as well but only to the extent o-f a quick 
qualitative look -for ccwnparison purposes. 

There is a great deal that can be done to enhance the use-ful- 
ness o-f the modeling equations. In particular, values for the 
K<P> function could be tabulated based on the examination of 
available propeller data for various combinations of governing 
parameters. For example, the first few terms of the Fourier 
series expansions could be tabulated for given pi tch- to-di a- 
meter and expanded area ratio combinations and used in a table 
look-up scheme such that a ship designer wishing to simulate a 
given propeller would need only select a particular rating 
point and combination of P/D and EAR ratios. The computer 
program could then select the necessary coefficients and pro- 
ceed with the computations. Or, if a comparative study of the 
rating points of various well designed propel 1 ers used in ship 
propulsion were to be conducted, a data base could be estab- 
lished such that a ship designer would need only select a 
generic type of propeller to begin a simulation. The computer 
would determine the appropriate rating point and correspondi ng 
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coe-f-ficients to initialize the program. 

The -form o-f the modeling equations make them particularly 
use-ful -for applications to controllable, reversible pitch 
propellers <CRPP) . Not only are they convenient -for conducting 
simulations, but they could -form the necessary -framework -for a 
control system synthesis. Such a control system could be 
use-ful in automated propulsion systems particularly involving 
CRP propellers. 

There are other potential uses including applications 
involving high thrust, low speed screws like those used in 
small harbor cra-ft or tugs. The applications to exotic 
propellers with very high skew rates or perhaps hybrid 
propellers are also worthy o-f investigation. Yet another use 
involves an examination of supercavi tati ng propellers and they 
modeling for simulation. Nith the promising results obtained 
here, follow— on efforts could be quite rewarding. 
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Appendix A 



Essentials o-f Eulerian Turbofnachines 



In a 'black-box' representation, a propeller like any turbo- 
machine can be modeled as in Figure Ai where the various 
parameters are torque < t) , rate o-f sha-ft rotation <N) , fluid 
velocity <V) , and -fluid -flow rate <F) . Ne de-fine an IDEAL 
propeller as one that, in the process o-f power conversion -from 
mechanical rotation to axial -fluid -flow, experiences no los- 
ses. Such a propeller may have little connection with reality, 
however the analysis o-f such a machine can provide qualitative 
insight into the mechanics o-f the energy trans-fer process. A 
more com-fortable model would include the e-f-fects o-f various 
types o-f losses into a PRACTICAL propeller model diagrammed as 
in Figure A2. 



The losses 


can be calculated by 


equation Al, 






Pd = T3^N 


- PXy 


<A1) 


where Pd 


represents the power 


di ssi pat i on . 


The amount o-f 



dissipation varies with the advance ratio such that two mini- 
mum will occur, one in the range o-f normal pump operation and 
the other approximately 180 degrees away in the -four quadrant 
diagram. Using the positive semi-de-f ini te dissipation -function 
described by Paynter C9] , the equation -for calculating the 
power dissipation -from the experimental data <in nondimension- 
al -form) is -found to be simply: 



&(P) = 2KCqVp - CtFp <A2) 
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TURBOMACHINE 






Figure A1 




Figure A2 
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The same results will be obtained using either a dimensional 
or a nondimensi onal formulation of the dissipation function. 
By plotting the square root of the dissipation versus the 
angle equal to the arctangent of J, the quadratic nature of 
the losses will result in a nearly linear graph such that the 
angle at which the minimum dissipation occurs is easily dis- 
cernabl e. 

Figure A3 shows the behavior of the dissipation as a function 
of the advance coefficient for propeller number 1. (Interest- 
ingly enough, of the two minimum values of 'del', the smaller 
occurs in the region of normal turbine operation.) The exten- 
sion to a PERFECT Eulerian turbomachine comes from proper 
scaling to eliminate the appropriate amount of dissipation so 
that the expression for the general dissipation function goes 
to zero at the prescribed operating point (best efficiency or 
minimum dissipation). The PRACTICAL Eulerian turbomachine then 
supplements this dissipation function by an amount equal to 
the minimum loss. 



Per Unit variables as used throughout this paper are defined 
as the nondimensi onal ratio of a parameter value to some 
predetermined base value. For a pump, those base values are 
taken to be at the point of best efficency or point of minimum 
dissipation. When presenting such data in a per unit form, 
obviously the reference operating point must be known to 
reconstruct the original experiment. 
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Figure A3 



Using per unit variables has significant advantages. First, 
they allow free manipulation of the applicable equations with- 
out regard to dimensionality. Secondly, the values of the per 
unit variables can be made to range in the neighborhood of 
easily handled numbers such as +1 and -1 by proper selection 
of the base values. Thus the per unit values give a qualita- 
tive as well as quantitative feel for their behavior. 

Since the propeller data which provided the basis for the 
derivation of the propeller equations was already in a non- 
dimensional form such that the physical dimensions of Minio- 
vich's propellers was no longer readily apparent, the analysis 
was done in a non-dimensional form to retain its generality. 
This may prompt one to conclude that the computational results 
are somewhat arbitrary, however one need only select Vp, N, 
and D to derive all the appropriate values. As utilized within 
the paper, the per unit values are found by: 

p=P/Po , t=T/To , n=N/No , q=F/Fo , etc. 
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Appendix B Selected Experimentally Determined Propeller 
Data Converted to ^ 4-Quadrant^ Form 

The data presented by Miniovich CIO] and converted to the 
modi-fied advance coe-f -f i c i en t -form Clll was converted to 4- 
quadrant -form by plotting Ct (thrust coe-f-f icient) and Cq 
(torque coefficient) in polar form utilizing the convention 
shown in Figure Bl. This is possible by recognizing that at a 
constant thrust or torque (which the Karman-Knapp Circle dia- 
gram utilizes), the defining relations used by Baker and 
Patterson shown below: 

Ct » T / pD*('^p2+N2D2) 

Cq = Q / PD3Cv^p2-»-N2D2) 

a = Arc tan •v'p/ND 
can be rewritten in per unit form as: 

R = 1 / SQRT( ABS(Ct) ) for the thrust curve 

and 

R = 1 / SQRT( ABS(Cq) ) for the torque curve 

and a = Arc tan q/n for both. 

As a result, the overlap of the data in the ranges of 
- 1 . 5< J<-0 . 667 and 0 . 667< 1/J< 1 . 5 act as a partial check by 
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plotting on top o-f each other. This kind oi representation of 
the data yields in two closed curves the complete characteris- 
tics of the propeller. Wylie and Streeter C223 address this 
method as applied to solution of transient problems using a 
digital computer. An additional advantage is the avoidance of 
the problems with singularities where ever propeller speed of 
advance or shaft rpm approaches zero. See also Reference C233. 
The curves are easily recognizible as typical for high speci- 
fic speed, axial flow pumps. See Reference C24] where Knapp 
also describes the 4-quadrant representation of pumps and 
turbines and outlines a procedure by which the representation 
could be used to deduce the transient behavior. Figure 62 is a 
sketch of the flow vs speed plane indicating the various 
regions of operation. The zero head and torque radial s are the 
asymptotes to which the thrust and torque curves approach for 
a given machine. 

Two series are shown here. Figures 63(a) and 63(b) are plots 
of thrust and torque for propellers 1,3, and 6. Figures 64(a) 
and 64(b) are thrust and torque curves for propellers 3,9 and 
15. The plots reveal the effect of variation in pitch-to- 
diameter ratio and expanded area ratio. 
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Figure BI 




Figure B2 
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Figure B3a 
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Figure B3b 
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Figure B4a 
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Figure B4b 
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Appendix C 



The Bond Graph Representation 



Paynter C9] presents the canonical bond graph representation 
•for an Eulerian turbomachine ass 



t 



n 



GY 

RG 



P 



q 



Rw 



Rq 



It is also represented as a modulated qyrator series adjoined 
to a shunt resistor: 

Rs 



r 



n 




In the -first case, the system equations are easily -found to 
be : 



m=R9Xq•^Rw3^n <D1) 



p = Rg X n - Rq X q <D2) 



using the pump sign convention. 
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In the second case, the equations are -found to be 



m = Rg 3^ q Rs <n-q) <D3> 

p = Rg 3C n - Rs X <n-q) <D4) 

By converting <D3) and <D4) to canonical -form, it is evident 
that Rs can be expressed as K<^> X v where 

V = SORT < nJ^n qJfq ) <D5) 



and 



Rg = R6 Rs X <n3fn q3fq)/v3«v 



<D6) 



A great deal o-f simplification is afforded with little loss in 
quality by making the simple approximation of Rg = n and Rs = 
K3^v where K absorbs the individual character of a particular 
machi ne. 

As a more widely applicable model, general dissipations could 
be adjoined to the 1-Junctions to account for fluid and mech- 
anical losses, i.e. 

Rs 
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The equations in matrix -form; 



'm‘ 




■Rs-*-Rdn 


Rg-Rs 




■ n ‘ 


P 

kh. J 




Rg-Rs 


Rs-Rdq 







<D7) 

<D8) 



A more complete model would leave the dissipation elements as 
•functions o-f the advance coe-f-f icient to absorb part o-f the 
individual behavior o-f a machine and smooth out the require- 
ments in the Rs term to -fit experimental data. 

The reader interested in pursuing the bond graph representa- 
tion is directed to Re-ference C25] -for an "entry level" ap- 
proach to bond graphs and to Re-ference C26] -for a more advan- 
ced treatment. 



73 



LIST of REFERENCES 



Cl] Kerwin,J.E., and Leopo1d,R., "A Design Theory -for Super- 
cavitating Propellers", Transactions SNAME, Vol . 72, 1964 

C2] Cummings,D.E. , “Numerical Prediction o-f Propeller Charac- 
teristics", Journal o-f Ship Research, SNAME Mar. 1973 

C3] Saunders,H.E. , "Hydrodynamics in Ship Design Mol. 1", 
Society o-f Naval Architects and Marine Engineers, 1957 

C4] Comstock ,J.P. , "Principles o-f Naval Architecture", Society 
o-f Naval Architects and Marine Engineers, 1967 

£5] Duncan, W.J., Thom,A.S. , and Young, A. D., "Mechanics o-f 
Fluids", Ptfnerican Elsevier Publishing Company, New York 1970 

C6] Tri tton ,D. J. , "Physical Fluid Dynamics", Mon Nostrand 
Reinhold Company, New York 1977 

C7] Csanady ,6.T. , "Theory o-f Turbomachines" , McGraw-Hill 
Series in Mechanical Engineering, 1964 

C8] Bovet,G.A., TRANS. ASME, 75; 975(1953) 

C9] Paynter ,H .M. , "The Dynamics and Control o-f Eulerian Turbo- 
machines", TRANS. ASME, Sept. 1972, pp. 198-205 

CIO] Miniovich ,Y.A. , "Investigation o-f Hydrodynamic Character- 
istics o-f Screw Propellers Under Conditions o-f Reversing and 
Calculation Methods -for Backing o-f Ships", BuShips Translation 
697, Translated by Royer and Roger, Inc. International 
Division, Nash. D.C. 1960 

Cll] Baker, D.N., and Patterson, C. L. , “Representation o-f Pro- 
peller Thrust and Torque Charac ter i st i cs -for Simulations", 
Naval Ship Research and Development Center, Annapolis Division 
Rpt MEL 202/67, Apr. 1968 

C 12] Rubis,C.J., and Moken,A.J., "Four Quadrant Open Nater 
Per-formance o-f Propellers in Terms o-f Modi-fied Advance Coe-f-fi- 
cients". Naval Ship Research and Development Center, Annapolis 
Division Rpt 6-175, May 1971 

C 13] Donsky,B. , "Complete Pump Characteristics and the E-f-fects 
o-f Speci-fic Speeds on Hydraulic Transients", Journal o-f Basic 
Engineering, TRANS. ASME, 1961 

C 14] Kovats,A. , "Design and Per-formance o-f Centri-fugal and 
Axial -Flow Pumps and Compressors", MacMillan Comp., New York 
1964 



74 



C 15] SI ocum,S.E. , "Speci-fic Speed Method of Propulsion Anal- 
ysis Naval Engineer's Journal, pg 727, Aug. 1954 

C 16] Rubis,C.J., “CRP Propel 1 er Ship Propulsion Dynamics 
Vol.l*, NSRDC rpt 6-101, 3238, Annapolis, Maryland Feb. 1971 

C 17] Peck,J.G., "Performance Characteristics of Three Prop- 
ellers Ni th ^^arying Pitch Distributions on an Inclined Shaft", 
NSRDC rpt SPD-497-02 Bethesda, Maryland Aug. 1974 

C 18] O'Br ien ,T.P. , "The Design of Marine Screw Propellers", 
London, Huthinson and Co. Ltd. 1962 

C 19] Wisl icenus,G.F. , TRANS. ASME, 78 : 1707< 1956) 

C20] MoKen,A.J., and Johnson ,R.D. , "Propulsion Dynamics Simu- 
lation of the SSN 685 Submarine Using a Di rec t-Curren t Elec- 
tric Drive", Naval Ship Research and Development Center, 
Annapolis Division, Rpt c27-823. Sept. 1974 

C21] Hecker,R., "Two-Quadrant Performance of Various Propel- 
l ers" , Naval Ship Research and Development Center, Bethesda, 
M.D., Rpt 386-H-35 Sept. 1973 

C22] Wylie,E.B., and Streeter , ,V.L. , "Fluid Transients", 
McGraw-Hill International 1978 

C23] MarchaljM., Flesh, G., and Suter,P., "The Calculation of 
Waterhammer Problems by Means of the Digital Computer", 
Proceedings, International Symposium Waterhammer Pumped 
Storage Projects, ASME, Chicago, Nov. 1965 

C24] Knapp, R.T., "Complete Characteristics of Centrifugal 
Pumps and Their Use in the Prediction of Transient Behavior", 
TRANS. ASME, Vol . 59, Nov. 1937, pp683-689 

C25] Karnopp,D., and Rosenberg , R . , "Analysis and Simulation of 
Multiport Systems, The Bond Graph Approach to Physical System 
Dynamics", M.I.T. Press, Cambridge, Mass. 1968 

C26] Paynter ,H.M. , "Analysis and Design of Engineering 
Systems", M.I.T. Press, Cambridge, Mass. 1961 



75 



1S94S7 



Thesis 

W58823 Wilken 
c.l An algebraic pro- 

peller model using 
Eulerian turbomachine 
theory . 



/ 



Thesis 199^^67 

W58S23 Wilken 

C'l An algebraic pro- 

peller model using 
Eulerian turbomachine 
theory. 



